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Abstract 

We introduce a speed selection mechanism for front propagation in reaction- 
diffusion systems with multiple fields. This mechanism applies to pulled 
and pushed fronts alike, and operates by restricting the fields to large 
finite intervals in the comoving frames of reference. The unique velocity 
for which the center of a monotonic solution for a particular field is insen- 
sitive to the location of the ends of the finite interval is the velocity that is 
physically selected for that field, making thus the solution approximately 
translation invariant. The fronts for the various fields may propagate at 
different speeds, all of them being determined though through this mech- 
anism. We present analytic results for the case of piecewise parabolic 
potentials, and numerical results for other cases. 



PACS numbers: 82.40. Ck, 05.45.-a, 47.54.+r. 



1 The selection mechanism 



In many systems rendered suddenly unstable, propagating fronts appear. The 
determination of the speed of a front propagating into an unstable state has at- 
tracted a lot of attention, since it cannot be achieved by simply solving an ordi- 
nary differential equation in the comoving frame of reference on a one- dimensional 
infinite domain. Indeed, there are many such solutions on such a domain, even 
though the propagating front in practice always relaxes to a unique shape and 
speed. The selection principles that have been formulated to determine the ob- 
servable front, without having to solve directly the partial differential equation 
of motion for a range of initial conditions, have involved concepts of linear and 
nonlinear marginal stability, of structural stability and of causality[l] , and all of 
them try to deal with the puzzle of the reduction of the multiple solutions to the 
unique observed one. All of these selection principles examine the waves from the 
viewpoint of the moving front, the corresponding wave equations being reduced 
then to ordinary differential equations involving the speed v of propagation. 

These various approaches can be problematic though in the case of multiple 
fields, because not all fields need propagate at the same speed, while the reduction 
of the set of partial differential equations to a system of ordinary differential 
equations requires that all fields be functions of the same variable x — vt. More 
recently, a complete analytical understanding of the propagation mechanism and 
relaxation behavior has emerged for those fronts that are "pulled along" by the 
spreading of linear perturbations about the unstable state, the so called "pulled" 
fronts [2] . This understanding resulted from a detailed study of the relevant partial 
differential equations and explains fully the behavior of pulled fronts. The speed 
selection mechanism for those fronts where linear analysis fails, the so called 
"pushed" fronts, is still however the subject of ongoing research for the case of 
multiple fields. 

The basic problem is that the ordinary differential equations that govern the 
motion of uniformly translating fronts do not seem to be able to determine the 



selected velocities for the various fields. Indeed, as we said above, the very exis- 
tence of different propagation speeds for the various fields makes the examination 
of the problem from the viewpoint of a particular moving frame of reference seem 
irrelevant. 



Let us take, for example, the equations 



dt dx 
d(p2 ^ 



01 - 01, 



D-^ + <t>2-<t>t + K<t>,, (1) 



dt dx^ 

where K is positive. The dynamics of 0i is always independent of that of 02, 
for fields propagating into the unstable state 0i = 02 = 0. II D <1, both fields 
move with speed v — 2. For > 1, the 0i and 02 fronts propagate with different 
speeds f i = 2 and V2 = '2\fD, respectively [2]. Clearly, the equations indicate 
that if 01 is a function oi x — vt then 02 should be too. Both fields should be 
propagating therefore with the same speed, which is the case however only for 
D <1. In fact, 02 always seems to be moving at the maximum available speed. 

It would appear thus that the fronts in reaction-diffusion systems with multi- 
ple components cannot be properly understood in terms of the properties of the 
ordinary differential equations that describe uniformly translating solutions. On 
the other hand, the examination of the full coupled partial differential equations 
is a rather complicated affair, and there is no universal way for dealing with 
pushed fronts. 

In this paper wc present a selection mechanism that applies to fronts invading 
both unstable and metastable states, whether they be pulled or pushed, and that 
works even for fields propagating at different speeds. Furthermore, it is easy to 
apply, since it involves examination of the system from the viewpoint of a single 
moving frame, resulting thus in coupled ordinary differential equations. 

This mechanism is the straightforward generahzation of the speed selection 
principle presented earlier for the case of a single field [3]. It exploits the fact that 
the observed front of a particular field is translationally invariant in the comoving 



frame of reference, even on a large finite interval, in the sense that its location 
is effectively independent of the ends of the interval. We shall be solving then 
the steady state equations of motion on a large finite interval with respect to a 
reference frame moving at an arbitrary given speed v, subject to the appropriate 
boundary conditions, obtaining a certain solution for each field. The solution 
for a particular field, however, will have approximate translational symmetry, 
thus becoming a physically observable front, only for a certain value of v. It is 
this value v* of v that is experimentally observed during the propagation of that 
field. Thus the selected front is the one that is ef fectively translationally 
invariant on a large finite interval, in the comoving frame of reference. Of 
course, this selected speed v* will not be appropriate, in general, for the other 
fields, in the sense that the corresponding solutions for the other fields need not 
have approximate translational symmetry at that speed. 

For values of v different from v* the midpoint of the field will be either at the 
left or the right end of the finite interval. It is only at v* that the midpoint can be 
anywhere near the center of the interval, becoming in fact indeterminate. Thus 
a graph of the speed v versus the midpoint of the field will have a plateau when 

V takes the value v*. The other fields will have such plateaus for other values of 
V. It is these plateaus, obtained from the graph of v versus the midpoints of the 
fields, that determine the physically selected speeds. 

Let us illustrate our mechanism with an example. We consider the following 
reaction-diffusion system: 

d4>i d'^(f)i dUi{(j)i) , , 

where each Ui is a function of the corresponding (f)i only. The fixed points of this 
system of differential equations provide the appropriate boundary conditions. 
Let us now assume that all fields have monotonic traveling wave solutions 
where — x — vt is the coordinate in a given moving frame of reference, with 

V > 0. Clearly, not all of these solutions need have translational symmetry. The 
above partial differential equations reduce then to the "steady state" ordinary 



differential equations 

d(j)i dUi{(f)i) 

We solve these equations on a large finite interval [Li,L2], with Li <g L2, subject 
to the boundary conditions (f>j{Li) — pj and 0j(I/2) — Qj, say, where j runs over 

all the fields, and where pj, qj and aij are constant. 

Let us now concentrate on a particular field (pi, and let us find the selected 
velocity of the corresponding front. Suppose that 4>i{^) is the solution of Eq. 
(3) subject to the boundary conditions mentioned above. There is only one such 
solution for a given velocity v. We multiply now Eq. (3) with d(j)i/d^ and integrate 
from Li to L2, obtaining thus 

with Wi{^) — d(j)i/d^. If is going to be a physically observable front on this 
large, but finite, interval, it will have to be essentially translationally invariant. 
This means that dcpi/d^ will be effectively zero in the regions close to the bound- 
aries, (pi having reached its fixed points there. Consequently Wi{Li) and Wi{L2) 
will tend to zero, while the integrals J^^ wf{^)d^ and jj^^ 'Wi{^)(j)j{^)d^ will be fi- 
nite and practically independent of Li and L2, as Li — > — 00 and L2 — > 00. Hence 
the speed v of Eq. (4) becomes independent of the endpoints of the interval, ac- 
quiring a unique value v*. In other words, only the front with that particular 
speed V* can correspond to an essentially translationally invariant field 0j. The 
other fields will not, in general, have translational invariance at that particular 
value of V, but that does not affect the above argument. Indeed, these other 
fields always appear multiplied by the quantity Wi{^), which is zero in the regions 
near the boundaries when v takes the value v* corresponding to a translationally 
invariant Thus the integral Jl^ '^iiO^'jiO^^^ remains independent of Li 

and L2, even if these other fields have no translational symmetry. 

The requirement that the front 0j be independent of the ends of the finite 



interval, when v — v*, selects therefore the speed 



(5) 



with Li — oo and L2 ^ oo, as the speed of the physically observed front for 
the field (pi. Note that no distinction has been made here between metastable 
and unstable states. Indeed, given any particular speed v, we can find a front 
interpolating between the stable and the unstable or metastable state, provided 
the solution is found on a finite interval. As the boundaries go to infinity, the 
value of the speed is restricted to v* and the front becomes the one corresponding 
to the speed of Eq. (5). 



2 Analytic example 

We shall demonstrate the proposed selection principle through analytic and nu- 
merical work. We shall adopt for our analytic work the following system of 
dimensionless partial differential equations: 







dt 








dt 




z/ > 3, and 






\u\ 




Ml- 



+ fu{<p2) + 9\<Pi\, (6) 



if \u\<l/2, 

I) if \u\>l/2. (7) 

This piecewise linear choice for the function /^(m) results from a piecewise parabolic 
potential and will lead to exact analytic solutions. Piecewise linear representa- 
tions of nonlinearities have often provided an analytically rigorous basis for the 
study of diffusion systems [4], as well as of nucleation and crystallization prob- 
lems [5], always on the interval (—00, 00). 

We shall be looking for uniformly translating solutions, functions of the vari- 
able ^ = X — vt, where > is an arbitrary given speed. Thus Eqs. (6) become 



d^(/) 



- + v^ + fM = 0, (8) 



and 



^^^ + v^ + M2)+g\cpi\ = ^- (9) 

The mirror symmetries of Eqs. (8) and (9) allow us to work with positive 
fields only. So we shall assume, without loss of generality, that 0i > 0, 02 > 0. 
When the fields are constant in space and time, they are at their fixed points. 
These fixed points (0i, ^2) are the points (0,0) (unstable fixed point), (0,1) (sad- 
dle point) and {l,l-\-g/v) (stable fixed point). We shall be interested in fronts 
invading the unstable state 0i = 0, 02 = 0, so we need to solve Eqs. (8) and (9) 
on the finite interval [Li, L2], subject to the boundary conditions 0i(Li) = 1, 
<f>2iLi) = 1 + g/u, 01 (L2) = 0, 02 (-^2) = 0, where Li <C <^ L2. These boundary 
conditions ensure that the system makes a phase transition from the unstable 
state to the stable state. Furthermore, we shall define the midpoints ^1 and ^2 
of the fields 0i and 02 through the relations 0i(^i) = 1/2 and 02(G) = 1/2, re- 
spectively, noting that the fields and their slopes have to be continuous at these 
points. 

The dynamics of 0i is decoupled from the dynamics of 02, consequently we 
can easily find the corresponding solution. There arc five boundary conditions 
that must be satisfied, namely two at the edges, the continuity of 0i and of d(j)i/d^ 
at ^1, and the definition of G. On the other hand, the solution of Eq. (8) for the 
field 01 will involve five unknown constants for any given value of namely ^ 
and two constants for each linear piece of /,i(0i). We expect therefore a unique 
solution 01 (^) for each value of v. 

Indeed, the exact solution of Eq. (8) for the field 0i is 



gmi(^-Li) _ gm2(^-ii) 

if ii<i< L2, (10) 



where 



k2^^{-v-^M^), (12) 



and ^1 satisfies 

The solution of Eq. (15) gives as a function of the speed v. We note that 
mi > > 1712 and k2 < ki < 0. For a given value of Eqs. (10)-(15) determine 
fully the field 01(0- 

As shown in our earlier work[3], the graph of v versus has a plateau at 
V = Vci, where Vd —{i-i + 1)/ \/2^ — 2. Indeed, if we require Li -C <^ L2, then 
Eq. (15) reduces to k2-\- mi — 0. This equation has a real solution, v — Vci, 
provided // > 3. In other words, for that particular value of v the midpoint ^1 
can be anjrwhere in the interval and cannot be determined, rendering thus the 
front effectively translationally invariant on the finite domain [Li, L2]. The value 
Vc\ is therefore the selected speed of 0i if > 3 (pushed case). If v > Vci, then 
^1 is close to Li, while for v < Vd we find that ^1 hes close to L2. The existence 
of the plateau at Vd shows thus that the only value of v for which the field 0i is 
approximately translationally invariant is Vd- 

Let us now turn our attention to the field 02- The relevant equation must 
be solved on the three segments determined by Li, L2, ^1 and ^2- There are 
seven boundary conditions: two at the edges, four from the continuity of 02 and 
d(j)2/d^ at the points ^1 and ^2, and the definition of {2- There arc also seven 
unknown quantities for a given value of v. two on each of the three segments, 
and ^2- Note that Eq. (15) determines the other midpoint. We expect therefore 
a unique solution 02(0 ^ach given value of v. 

Clearly, if g were zero, then the speed of propagation for 02 would be Vc2 — 
{u + l)/\/2u — 2, and the two fields would be completely decoupled. We need 
to examine what happens for g ^ 0. We can obtain some qualitative results by 



looking at Eqs. (8) and (9). We distinguish two cases, 
(a) C£ise Vc2 > Vci 

(i) We examine the case v > Vc2 > Vd first. Since v > Vd, we have ~ Li, 
as discussed in our earher work [3]. Therefore the field 4>i falls very rapidly from 

1 to 0, and it remains equal to on almost all of the interval [Li, L2]. Since 0i 
is approximately zero almost everywhere, Eqs. (8) and (9) decouple. Thus 02 
behaves as if g were equal to 0. Since v > Vc2, this implies that ^2 ~ -t'l- Hence (f>2 
is also zero practically everywhere. Both fields are essentially on the fixed point 
(0,0). 

(ii) We examine the case Vc2 > v > Vd next. Again ^1 ^ Li, since v > Vd, 
hence (f)i is essentially zero on almost all of the interval. The fields decouple once 
more, but now ^2 ~ -^2, since v < Vc2- Hence the field 02 is nonzero up to the 
point ^ — L2. The only fixed point that is available for the two fields is then the 
point (0,1). Thus 0i starts out at Li having the value 1 and very rapidly drops 
down to zero. The field 02, on the other hand, starts out having the value l+g/u 
at Li, drops down to the value 1 almost immediately, and then it stays there till 
it reaches the other edge, where it drops down to zero. Thus the midpoint of 

02 shifts abruptly from the left edge to the right edge the very moment we pass 
from case (i) to case (ii), i.e. at ^; = Vc2, because ^2 ~ -^1 for v > Vc2, but ^2 ~ -^2 
for V < Vc2- 

(iii) We examine finally the case Vc2 > Vd > v. Since Vd > v, we have ^1 ~ L2. 
Therefore 0i remains on the value 1 on almost all of the interval, dropping down 
to zero only very close to the right edge. The only available fixed point for the 
two fields is then the point (1, 1 + g/i^)- Therefore the field 02 remains stuck 
at the value g/v almost everywhere, dropping to zero only very close to the 
right edge, whereby ^2 ~ -^2- We note that the abrupt shift of the midpoint of 
01 occurs ak, V — Vd- 

These arguments indicate then that the field 0i acquires approximate trans- 
lational invariance when v — Vd, as expected, since its dynamics is decoupled 



from the dynamics of 02- At that speed we have a plateau of v versus ^i. On 
the other hand, the plateau of v versus ^2 occurs ai v = Vc2- Therefore the front 
of 02 propagates with the speed Vc2- Indeed then, our finite interval mechanism 
gives the selected velocities for both fields, in spite of their being different. We 
also note that in all likelihood ^1 < ^2, since in the region Vd < v < Vc2 we found 
that ^1 ^ Li and ^2 ~ L2. 

(b) Ceise Vc2 < Vd 

(i) We examine the case v > > Vc2 first. Since v > Vd, we have ^1 ~ Li, 
as discussed in our earlier work [3]. Therefore the field 0i falls very rapidly from 

1 to 0, and it remains equal to on almost all of the interval [Li, L2]. Since 0i 
is approximately zero almost everywhere, Eqs. (8) and (9) decouple. Thus 02 
behaves as if g were equal to 0. Since v > Vc2, this implies that ^2 ~ -^i- Hence 02 
is also zero practically everywhere. Both fields are essentially on the fixed point 
(0,0). 

(ii) We examine the case Vd > v > Vc2 next. Now ^1 ~ L2, since v < Vd, 
hence 0i is stuck on the value 1 on almost all of the interval. The only fixed point 
that is available for the two fields then is the point {1,1 + g/u). That means that 

02 must be stuck at the value l + g/v on almost all of the interval, dropping down 
to zero only close to the right edge. Hence ^2 ~ -^2- Thus both the midpoints of 
02 and 01 shift suddenly from the left edge to the right edge the very moment 
we pass from case (i) to case (ii), i.e. at v = Vd- 

(iii) We examine finally the case Vd > Vc2 > v. Since Vd > v, we have ^1 ~ L2. 
Therefore 0i remains on the value 1 on almost all of the interval, dropping down 
to zero only very close to the right edge. The only available fixed point for the 
two fields is once more the point {1,1 + g/u). Therefore the field 02 remains stuck 
at the value 1 + g/u almost everywhere, dropping to zero only very close to the 
right edge, whereby ^2 ~ -^2- 

These arguments indicate then that the field 0i acquires approximate transla- 
tional invariance when v — Vd, since its dynamics is decoupled from the dynamics 



of 02- At that speed we have a plateau of v versus ^i. On the other hand, the 
plateau of v versus ,^2 occurs also at v = Vd- Therefore the front of 02 propagates 
with the speed Vd- In this case both fields propagate at the same speed. 

We note that the field 02 always propagates at the maximum possible speed, 
i.e. Vc2 in case (a) and Vd in case (b), just like the fields of Eqs. (1). Our finite 
interval mechanism is able to handle both cases though. We also note that the 
speeds of propagation are independent of the coupling constant g, irrespective of 
how large or small it is. 

Let us now verify this behavior by solving analytically Eq. (9) to find 02, 
given the solution of Eq. (10) for the field 0i. 

We shall assume that ^1 < ^2, for the sake of definiteness. This situation is 
appropriate for the case Vc2 > Vd, according to the arguments presented above. 
If contradictions arise due to this assumption, it will be easy enough to repeat 
the work with the contrary assumption. In fact, it turns out that the relation 
Ci ^ holds even if Vc2 < Vd, in the examples we shall present. 

Let us examine the region ^2 < C ^ -^2 first. The boundary conditions 
02(S) — 1/2 and 02(-^2) = determine the two constants that will appear in 
the solution of the ordinary differential equation on this interval. Thus the full 
solution for 02 on the interval [^2, L2] turns out to be 



where 



(17) 



and 



i - (6 - L2)(-e'=i«2-^2)/(2A;i + v) + e'^^^^^-'^^y {2k2 + v))zi 



(18) 




^ 2ki + v 2k2 + v ^'^ e-'=2(€2-i:'2) _ e-fci(«2 

UQk-iii2-L2) _ ^ gfc2(6-i2) 



(19) 



We can now use the known values of 02(S) ^-iid ^2(^2) as boundary conditions in 
order to solve Eq. (9) on the interval [^1,^2]- We find 

Me) = 1 + ^ + Ae"i(^-^2) + Se"2(^-^^\ (20) 

u + 1 



where 



ni = ^{-v + Vv^ + Au), (21) 
^2 = ^(-^ - Vv^ + 4:iy), (22) 



with ni > and n2 < 0, 



^ 1 _ _^M2-L2) _ (23) 

^2 = 0;(6) - ^(^le'^^^^-^^) - fc2e'=^(«^-^^)), (24) 



-ni ($2-^2) ^2 - ^2!^! 



Til - 77,2 
^2 - Til 

We can now use Eq. (20) to find ^2(^1) ^.nd 02(^i)- We get 



(25) 
(26) 



(/,2(A) = 1 + + Ae"i(^i-^2) + 5e"2(^i-^^^ (27) 

and 

</,;(^^) = (A;i 6^=1 (^1-^2) - A;2e'=2(^i-^2))-^ +niAe"i(«i-^2) +n2Se"2(^i-^2). (28) 

Finally, we can use these values of 02 and 02 at ^1 in order to find the solution to 
02 in the interval [Li, ,^1]. Imposing then the boundary condition 02(i^i) = l+g/u 
on this solution leads to the final relation 

Q2 — TliQi = 

+ (Q2 - n2Jli)e("i-"2)(-^2+-^i) 

+ f- + ^ - TT^) (n2e-("i-"^)(^i-^i) - ni)e"2(«^-^i). (29) 



We can use this equation to find the plateau of 02- 

Indeed, let us take the case Vc2 > Vd first. If v > Vd, then we must have 
^1 ~ Li. In that case we can show that any ^2 that is far from Li and L2 will 
satisfy Eq. (29), provided A;2 + rii = 0, a relation equivalent to the requirement 

that V be equal to Vc2- Thus 02 has a plateau above Vd, at the speed v = Vc2- 
Note also that the midpoint ^2 of 02 is already at L2 when the midpoint of 0i 
shifts to L2. 

On the other hand, if we look at the case Vc2 < Vd, then we see that below 
Vd we must have ^1 ~ L2. But since all the above analytic equations have 
been derived under the assumption ^1 < ^2, we conclude that ^2 must be close 
to L2 as well. Thus, when the midpoint ^1 of 0i shifts abruptly to the right 
edge, it forces the midpoint of the other field to go there as well, provided the 
analytic equations have solutions consistent with the assumption ,^1 < ,^2- The 
corresponding behavior is illustrated by the examples of Figs. 1 and 2, verifying 
thus the qualitative conclusions drawn earlier. In particular, these figures confirm 
the assumption ^1 < ^2- 

We see then that, if Vc2 > Vd, the field 02 has a plateau at the highest speed 
Vc2- If, on the other hand, Vc2 < Vd, then 0i pulls 02 and forces it to propagate 
at the higher speed Vd- This behavior is seen even when the coupling constant 
g takes very small values, and matches the behavior of the fields that obey Eq. 



3 Numerical examples 

We can demonstrate our selection mechanism numerically as well. Let us examine 
the following system: 



(!)• 



901 



5^01 



dt 

902 



^202 



+ hb^{(j)2) + gh, 



(30) 



dt 



where 

hb{u) = ^u{l-u){b + u). (31) 

It was this particular choice of h{u) that was used for the case of a single field when 
the concepts of linear and nonlinear marginal stability were first proposed [6]. 
That study found that for < 6 < 1/2 the selected speed of the single field for 
the front invading the unstable state is (26 + l)/\/2b. 

We shall consider values of b less than 1/2 (pushed case). Thus, if the cou- 
pling constant g were zero, then the two fields (pi and (p2 of Eq. (30) would 
propagate separately, with speeds = (26i + 1)/ \/26i and Vc2 = (262 + 1)/ \/262, 
respectively. 

We have solved Eqs. (30) numerically on a finite ^ domain for the h{u) of 
Eq. (31), assuming that both fields are functions of the variable ^ = x — vt, with 
g=21, bi — 1/8 and 62 — 1/18, subject to the boundary conditions (f)i{Li) — 1, 
0i(-^2) = 0, 02(-^2) = and 02(-^i) = 1-5. These values correspond to the stable 
and unstable fixed points (1,1.5) and (0,0) for the fields 0i and 02, the points 
in other words where Eqs. (30) acquire uniform solutions. The solutions that 
interpolate between these two fixed points are shown in Fig. 3. We can see that 
the two fields are at the saddle point 0i — 0, (j)2 — 1, on almost all of the interval. 
This feature reminds us of the dual fronts, where the decomposition from the 
unstable to the stable state proceeds via an intermediate saddle point, with the 
two fields propagating at different speeds [7]. 

Fig. 4 shows the locations ^1 and ^2 where the fields 0i and 02 attain the 
value 1/2, at a given arbitrary speed v, when g = 21, bi = 1/8 and 62 = 1/18. 
We can deduce from this figure that the field 0i propagates with speed v^i = 2.5, 
while the field 02 propagates with speed Vc2 — 10/3. Indeed, we see that at these 
velocities there are plateaus of v versus the location where each field acquires 
the value 1/2. In other words, the location of the midpoint of the front for the 
field 01 or 02 is indeterminate at the corresponding velocity Vd or Vc2, rendering 
the solution essentially translation invariant there. We note that the two fronts 



propagate at different speeds, and that < ^2- 

We have also solved Eqs. (30) for the case = 21, 61 = 1/18 and 62 — 1/8. 
The unstable fixed point is still the point (f)i — (p2 — 0, but the stable fixed point 
is the point 0i — 1, (f)2 — 1.7768, since these values satisfy Eqs. (30). Hence the 
solutions we seek have to interpolate between these two fixed points, subject to the 
boundary conditions 0i(i^i) = 1, 0i(-Z^2) = 0, 02(-^2) = and (f>2{Li) = 1.7768. 
Fig. 5 shows the midpoints ^1 and ^2 of the two fields at an arbitrary speed v. 
We observe once more that there are the usual plateaus, ^1 being again less than 
^2. However, both plateaus occur at the speed v — — 10/3. The field 02 
is pulled by (f)i and is forced to propagate at v^i, rather than at its own lower 
speed Vc2 = 2.5. The existence of the common plateau indicates once again that 
at that particular speed Vd the locations of the midpoints of the fields become 
indeterminate, making thus the two fields effectively translation invariant. 

4 Concluding remarks 

We see then that requiring the solution of a field to have approximate translational 
invariance on a finite interval in the comoving frame of reference results in the 
selection of a speed for the front. We can adopt then a selection principle that 
reads "the selected front is the one that is approximately translationally invariant 
on a large finite interval, with respect to the comoving frame of reference" . This 
principle is very easy to implement, especially numerically. Indeed, it suffices 
to solve the moving frame equation on a large finite interval, for an arbitrary 
propagation speed. For large speeds we expect the midpoint of the front to be 
close to the left boundary. As the speed is lowered, the midpoint suddenly moves 
to the right boundary. The speed v* at which this sudden move occurs is the 
speed selected by the physically observed front. 

Of course, the selection of a certain plateau for a particular field need not 
include the selection of such a plateau for other fields as well. In fact, the advan- 
tage of our mechanism is that the graph of the speed of propagation versus the 



midpoint of each field allows us to find the plateaus for all the fields, and hence 
the corresponding physically selected velocities. The approximate translational 
invariance that each field acquires as its plateau is reached need not involve the 
other fields as well. 

Thus our mechanism can find the selected velocities for all the fields present 
in a system of coupled partial differential equations, by simply solving a much 
simpler system of coupled ordinary differential equations, having assumed that 
all the fields are travehng at the same arbitrary speed v. The value of v for which 
the midpoint of a particular field becomes indeterminate is the physically selected 
velocity for that field. 
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Figure Captions 



Figure 1: The speed as a function of the midpoints of the fronts invading the unsta- 
ble state, for the system of Eqs. (8) and (9), with Li — —10 and L2 — 10, 

g = 9.5, /i = 9, z/ = 19. The plateau is at f = 2.5 for 0i, and at f = 10/3 
for 02- AH quantities are dimensionless. 

Figure 2: The speed u as a function of the midpoints of the fronts invading the unsta- 
ble state, for the system of Eqs. (8) and (9), with Li = —10 and L2 = 10, 
g — 4.5, fjL — 19, u — 9. The plateau is at = 10/3 for both 0i and 02- All 
quantities are dimensionless. 

Figure 3: The solutions of Eqs. (30) that interpolate between the fixed points (0,0) 
and (1,1.5), for Li = -15, L2 = 15, -u = 2.7, g = 21, bi = 1/8, 62 = 1/18. 
These solutions he on the saddle point (0,1) on most of the interval. 

Figure 4: The speed v as a function of the midpoints of the fronts invading the unsta- 
ble state, for the system of Eqs. (30), with Li — —15 and L2 — 15, g — 21, 
61 = 1/8, 62 = 1/18. The plateau is at = 2.5 for 0i and at v — 10/3 for 
4>2- All quantities are dimensionless. 

Figure 5: The speed v as a function of the midpoints of the fronts invading the unsta- 
ble state, for the system of Eqs. (30), with Li = —15 and L2 = 15, g = 21, 
bi = 1/18, 62 = 1/8. The plateau is at v = 10/3 for both 4>i and 4>2- All 
quantities are dimensionless. 



